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摘 S. 在 太阳 观测 研究 中 ， 高 分 辨 图 像 与 全 日 面 像 的 配 准 是 一 项 非常 有 意义 的 工作 
但 由 于 他 们 之 间 存 在 着 旋转 、 缩 放 和 平移 ,因此 很 难 高 精度 地 进行 匹配 。 提出 一 种 结合 局 部 
统计 信息 和 控制 点 匹配 的 图 像 配 准 方法 , 核心 思想 为 将 视 场 等 间隔 划分 为 大 量 重要 的 局 部 区 


域 , 通过 相关 匹配 在 全 日 面 像 上 寻找 对 应 的 局 部 区 域 , 然后 计算 每 一 对 局 部 区 域 间 的 亚 像 元 


据 控 制 点 建立 仿 射 变换 的 转换 方程 , 采用 最 小 二 乘 求解 整个 视 场 的 转换 参数 , 根据 解 出 的 参 
数 重 新 对 图 像 进行 和 迭代， 得 到 收敛 后 的 结果 并 进行 配 准 。 通 过 对 高 分 辨 观测 图 像 和 全 日 面 
SDO/HMI 连续 谱 图 像 进行 配 准 ， 拟 合 结果 的 偏差 在 0.25” 以 内 。 


关键 词 : 太阳 图 像 配 准 ; ”高 分 辨 图 像 ; 全 日 面 图 像 ; 
中 图 分 类 号 : P182.2 ”文献 标识 码 : A 文章 编号 : 1672-7673 


在 研究 太阳 图 像 的 观测 数据 时 , 经 常 需要 将 不 同 设备 的 观测 数据 进行 对 比 ， 从 而 发 挥 它 
们 各 自在 不 同 波段 、 不 同 视 场 、 不 同 分 辩 率 上 的 观测 优势 。 太 阳光 学 成 像 观 测 包 括 全 日 面 和 
高 分 辨 观测 两 类 。 对 于 全 日 面 图 像 ， 由 于 太阳 几何 尺寸 和 运动 方式 已 知 ， 如 太阳 的 半径 、P0 
角 、 自 转 等 ， 因 此 各 种 观测 设备 均 把 自己 的 数据 在 预 处 理 后 归 算 到 标准 的 太阳 坐标 系 。 然 而 
对 于 局 部 高 分 辩 像 ， 由 于 观测 设备 的 指向 误差 、 视 场 旋转 以 及 焦距 的 变化 ,使 得 直接 依赖 仪 
器 参数 确定 准确 的 太阳 视 场 坐标 是 不 精确 的 , 因此 只 能 借助 太阳 图 像 上 相似 的 太阳 结构 (如 
光 球 层 对 应 光 球 层 ， 色 球 层 对 应 色 球 层 ) 与 全 像 进行 转换 ， 即 视 场 匹配 。 然 而 全 日 面 观 
测 有 大 于 31” 的 视 场 , 但 最 高 的 空间 像 元 分 辨 仅 为 0.5”《 和 空间 太阳 动力 天 文 台 的 日 球 层 磁 
场 观测 仪 ，SDO/HMI)。 而 地 面 高 分 辨 观测 ， 如 中 国 的 1m 新 真空 太阳 望远镜 (New Vaccum 
Sloar Telescope,NVST) 和 美国 大 能 湖 天 文 台 的 新 太阳 望远镜 (New Solar Telescope,NST), 视 
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场 仅仅 为 1 ”到 2”， 但 像 元 分 辨 率 达 到 0.05” 以 上 ， 这 就 意味 着 需要 将 两 个 空间 分 辨 率 相 
差 十 倍 的 图 像 进行 配 准 。 即 便 不 考虑 其 他 畸变 ， 局 部 像 的 视 场 旋转 R EARN (S) 和 中 
Ui Bb CD Ty) 是 必须 得 到 的 4 个 参数 。 
图 像 配 准 是 指 在 不 同时 间 、 不 同 条 件 下 得 到 的 并 且 位 于 不 同 坐 标 系 下 的 同一 场景 的 两 幅 
或 者 多 幅 图 像 进行 匹配 的 过 程 "。 由 于 存在 各 种 不 同 的 情况 ,同一 场景 的 多 幅 图 像 会 在 分 辨 
率 、 位 置 、 尺 度 、 成 像 模式 等 方面 存在 很 多 差异 ， 图 像 配 准 是 找到 一 个 几何 变换 的 最 优 解 ， 
使 得 图 像 上 相对 应 的 点 , 或 者 至 少 是 有 意义 的 点 达到 匹配 状态 。 目 前 在 图 像 配 准 领域 , 常见 
的 配 准 方法 有 基于 特征 的 点 匹配 方法 ， n 图 像 中 的 特征 点 , 或 者 显著 特征 构造 描述 子 确 
定 的 匹配 点 , 然后 构造 转换 方程 进行 配 准 。 理论 上 这 个 方法 如 果 有 足够 多 和 足够 精确 的 控制 
点 是 可 以 很 好 地 解决 两 个 图 像 的 几何 变换 ， em 缩放 、 平 移 等 问题 。 但 关键 问题 是 如 何 
确定 这 些 点 ， 这 方面 在 图 像 处 理 界 有 很 多 研究 ， 如 Harris 角 点 及 其 衍生 的 Harris-Laplace 77 
法 外 ， 尺 度 不 变 特征 (Scale-invariant feature transform,SIFT) 号 及 其 衍生 的 PCA-SIFT 算法 
(Principal Com-ponent Analysis SIFT)“, au P(E FBI. RRR EAB PIE BE (Very Fast 
SIFT, VF-SIFT) ", WESE (Speeded Up Robust Features, SURF) “等 。 另 外 一 
种 是 基于 区 域 信 ) 息 的 统计 配 准 方法 , 主要 通过 图 像 区 域 的 统计 信息 度量 图 像 的 相似 程度 ， 如 
互信 息 ”， 互 相关 函数 "”"”， 模 板 匹 配 “"， 不 变 矩 ， 相 关 匹 配 *"“ 等 。 但 一 般 来 讲 这 种 方法 
主要 解决 图 像 间 的 平移 问题 对 于 具有 旋转 和 比例 缩放 的 图 像 ， 又 引申 出 在 频率 域 的 相位 相 
eR, Fouier-Mellin 5 |“ , 
然而 ， 由 于 太阳 观测 图 像 和 普通 自然 图 像 有 很 大 的 区 别 ， 其 中 的 特征 边界 不 清晰 ， 结 构 
相对 模糊 且 存 在 很 多 噪声 , 同时 动态 范围 大 , 运动 速度 较 快 , 并且 存在 很 多 不 同 的 活动 现象 ， 
大 部 分 运动 现象 为 非 刚性 运动 ; 同时 太阳 图 像 的 图 幅 一 般 较 大 ,有 效 目 标 小 且 缺 少 有 效 的 纹 
H, 不 存在 明显 的 梯度 变化 等 特征 , 并 且 由 于 拍摄 太阳 图 像 的 设备 以 及 拍摄 过 程 中 的 不 确定 
性 ， 采 集 的 太阳 图 像 还 存在 旋转 、 缩 放 及 平移 等 情况 ,不 同 图 像 的 空间 分 辨 率 经 常 不 在 同一 
尺度 ， 这 就 为 确定 特征 对 应 点 带 来 非常 大 的 困难 。 长 期 以 来 , 在 同一 观测 设备 的 太阳 时 间 序 
列 图 像 的 视 场 对 齐 一 直 采 用 基于 互相 关 的 统计 配 准 ， 而 不 同 观测 设备 之 间 的 图 像 配 准 , 还 停 
留 在 利用 一 些 显著 的 太阳 特征 ， 如 黑子 、 气 孔 、 大 米粒 构造 对 应 点 的 方式 。 然 而 这 些 特 征 点 
数量 少 , 结构 不 规则 , 难以 保证 位 置 测量 精度 , 同时 也 难以 使 用 计算 机 技术 进行 自动 化 测量 。 
尤其 是 将 高 分 辨 局 部 像 对 应 到 全 日 面 像 的 时 候 ， 这 种 问题 更 为 突出 。 但 从 另 一 方面 ,通过 望 
远 镜 机 械 和 终端 光学 系统 ， 对 局 部 太阳 像 指 向 、 比 例 尺 和 位 置 有 初步 估计 ， 只 是 这 些 数 值 不 
e 很 精确 且 随 时 间 有 一 些 变化 。 因此 实际 面临 的 是 如 何 提 高 配 准 精度 ,进行 - 图像 匹配 的 问 


题 。 


3g 


m 


基于 这 些 问题 ， 本 文 提 出 了 一 种 将 统计 信息 和 相关 点 匹配 相 结合 的 太阳 图 像 配 准 算 法 ， 
解决 将 不 同 分 辨 率 下 的 局 部 高 分 辨 观测 数据 和 全 日 面 数 据 进 行 高 精度 坐标 转换 的 问题 , 其 中 
包括 测量 图 像 指向 角 、 像 元 比例 尺 和 图 像 中 心 在 全 日 面 像 上 的 位 置 。 这 一 方法 的 基本 思想 是 
在 初步 匹配 的 基础 上 , 首先 将 视 场 划分 成 大 量 重生 的 局 部 区 域 代 蔡 点 匹配 中 的 控制 点 , 通过 
基于 统计 的 相关 匹配 确定 对 应 区 域 间 的 亚 像 元 位 移 , 从 而 确定 特征 点 对 在 同一 坐标 系 中 的 位 
Bu 然后 用 最 小 二 乘 欠 代 求解 整个 视 场 的 转换 参数 ， 并 逐步 迭代 提高 拟 合 的 精度 ， 成 功 将 
1M 太阳 望远镜 和 新 太阳 望远镜 的 Ti0 观测 数据 与 几乎 同时 刻 《〈 最 长 间隔 11 s) 的 SDO/HMI 
连续 谱 进 行 了 匹配 ， 拟 合 结果 的 偏差 可 达到 0. 25” 之 内 。 


1 实验 数据 来 源 


本 文 实验 使 用 的 全 日 面 图 像 来 自 太 阳 动 力学 观测 卫星 SDOCSolar Dynamics Observatory, 
SDO) 上 日 球 层 磁 场 观 测 仪 HMI (Helioseismic and Magnetic Imager, HMI) 采集 的 连续 谱 全 
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日 面 像 ， 而 高 分 辨 图 像 为 1M 太阳 望远镜 和 新 太阳 望远镜 的 TiO FERRER. SDO/HMI ik 
的 空间 分 辨 率 约 为 0. 5043 arcsec/pixel， 图 像 尺 寸 为 4096 像 元 *4096 Rac; IM 太阳 望 远 
镜 和 新 太阳 望远镜 图 像 的 像 元 比例 尺 为 0. 04 arcsec/pixel 和 0. 03 arcsec/pixel。 

图 1 (a) 为 IM 太阳 望远镜 的 Ti0 观测 数据 ， 观 测 时 间 为 2013-7-15 07:29:09， 有 效 
区 域 尺寸 为 2304 像 元 *1920 像 元 ; 图 1 (b) 为 SDO/HMI 连续 谱 观 测 数据 ， 观 测 时 间 为 
2013-7-15 07:29:15, 红 框 标记 区 域 为 IM 太阳 望远镜 在 SDO/HMI 中 的 大 体位 置 。 本 文 主 
要 以 此 图 像 的 匹配 说 明 算法 流程 。 


a IM 太阳 望远镜 (Tio) b SDO/HMI 
图 1 观测 数据 


Fig. 1 Observed data 


2 配 准 的 方法 


2.1 配 准 方法 的 流程 
配 准 方法 的 主要 流程 如 图 2。 
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开始 


估计 RST 参 数 
Y E aR 
旋转 SD0/HMI 降 采 样 NVST 像 
-— — 
高 斯 平滑 
v v 
截取 子 图 1, 划分 N 个 子 块 大 


计算 六 和 及 互相 关 


测量 六 和 fi 亚 像 元 偏 移 


y 
ETCHEA Pi Pr 


Y 


建立 转换 方程 


A A 
最 小 二 乘 求 解 RST 参 数 


未 收敛 


收敛 
图 2 算法 流程 


Fig.2 Algorithm flow chart 


图 2 可 以 描述 为 下 列 11 个 步 又: 
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(1) 对 IM 太阳 望远镜 像 的 指向 、 比 例 尺 和 位 置 进行 初步 估计 ， 得 到 旋转 〈Rotation)、 缩 


W (Scaling) KFE (Translation) 参数 (下 文 称 RST 参数 ); 


(2) 根据 RST 参数 ， 对 SDO/HMI 全 日 面 像 进行 旋转 ， 并 截取 与 1M AP 
的 SDO/HMI F BII, ; 


望远镜 像 视 场 相 同 


(3) 根据 RST 参数 和 SDO/HMI 的 比例 尺 , 将 1M 太阳 望远镜 像 缩小 至 和 SDO/HMI 同一 比例 尺 


下 ， 即 对 1M 太阳 望远镜 图 像 进 行 降 采样 ; 


(4) 对 降 采样 后 的 1M 太阳 望远镜 像 进行 高 斯 平滑 , 减 小 采样 过 程 中 的 频率 混合 。 平滑 后 的 


IM 太阳 望远镜 降 采 样 像 用 Ik 表示 ; 


(5) FENZ I 等 间隔 划分 为 YX 个 重 琶 的 子 块 ， 并 将 每 个 子 块 记 为 六 ; 


(6) ”对 每 一 个 子 块 六 和 SDO/HMI. 子 图 /利用 归 一 化 互相 关 方 法 求 得 大 对 应 的 每 一 个 


TY f 


C “对 每 一 组 大 和 六 ， 测 量 其 亚 像 元 偏 移 ， 并 将 子 块 中 心 设置 为 一 组 对 应 的 控制 点 ; 
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(8) 筛选 每 组 控制 点 ， 将 控制 点 间距 离 大 于 一 定 像 元 的 点 剔除 ; 
(9) 根据 剩 下 的 控制 点 ， 建 立 仿 射 变换 的 转换 方程 ， 采 用 最 小 二 乘 求解 RST 参数 ; 


C100 计算 残 差 ”和 选 代 次 数 了 上 ， 若 满足 { <05^ <01}; [ <03}; 


{L>30, =min( ，aj 任 一 条 件 ， 则 停止 迭代 ， 和 否则 ， 返 回 步 又 (2)， 根 据 步骤 CoD 解 得 


的 RST 参数 重新 对 SD0/HMI 和 1M 太阳 望远镜 像 进 行 迭 代 处 理 ; 
(11) 最 终 得 到 迭代 完成 后 的 RST 参数 并 进行 配 准 。 
在 上 述 步 骤 中 ， (60 至 COD 步 为 算法 的 核心 步骤 ,下 文 对 该 部 分 进行 更 为 具体 描述 。 
需要 注意 ， 开始 适 代 后 (2) 到 CO 步 中 的 RST 参数 会 根据 上 次 大 代 的 结果 进行 下 一 次 送 
代 ， 意 味 着 算法 会 重新 对 图 像 进 行 平滑 
2. 2 分 块 确定 特征 控制 点 


WIS 为 待 匹配 图 像 ,为 参考 图 像 。 将 Lr 划分 N 个 特征 区 域 , WIRE ft 912. N) 


e I BURGTZSU V , SET DOR fL BINE 2gu v ， 可知 u<U,v<V ， 则 参考 图 像 1 和 图 像 fr 


HI EHK K% Corr (x, y) 定义 为 ， 


SG») FA 7] 
Corr(x,y) = ——— TTE = (D 
Po (xs) 大 | $6) f 
= Joh, (st) nA Liebe, KeA]: n-En] G9 


1, All fi EAB IRA AP ELA, RIER, 和 f] TRES (9,1) 进行 。 


根据 上 述 公 式 ， 每 个 大 和 到 做 归 一 化 互相 关 ， 可 得 到 一 个 互相 关 曲 面 ， 而 峰值 的 坐标 


o 代表 了 fi 图 像 中 心 在 1 中 的 整数 像 元 精度 的 匹配 位 置 。 


在 I, 中 根据 上 述 位置 可 截取 和 同样 大 小 的 特征 区 域 ,得 到 每 一 个 大 对 应 的 fi. BI 


们 组 成 的 集合 记 为 {fi,i=1,2...,N}。 


对 于 两 组 对 应 的 特征 区 域 的 fe M fso EMRE F 描述 ff. 任意 两 个 特征 区 域 之 间 
的 关系 ， 即 


Rap(A RIN (2) 


F, 表示 任意 两 个 对 应 的 局 部 特征 区 域 ( 太 , Fi) 。 


对 于 每 一 对 特征 区 域 ( 友 , 太 ), 计算 标准 的 互相 关子 数 。 假设 DD 为 特征 区 域 ( 友 , 太 ) 间 
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HB B c^) Fre LR RE DC, PR ETE RET S P EG Ap (m, n) , 


Sex) XQUY) 
xut xn 


dich, [AHR WU I' I —b, b=min(7),1eD 。 tiz (fs. fs) 图 像 间 互相 关 函 数 的 


， (3) 


局 部 重心 坐标 代表 每 一 对 特征 区 域 ( 大 , f) 间 的 亚 像 元 偏 移 。 


根据 亚 像 元 偏 移 的 情况 ， 可 以 建立 每 一 对 (万 , 友 ) 图 像 中 心 点 对 应 的 坐标 , 得 到 和 大 之 
间 的 对 坐标 转换 控制 点 。 用 集合 P 描述 控制 点 对 之 间 的 关系 ， 即 


P=P(pi.ps).1Si<N . (4) 
2.3 最 小 二 乘 求解 转换 参数 
考虑 平移 、 旋 转 和 缩放 的 情况 ， 假 设 控 制 点 对 (pi, p;) ， 对 应 的 坐标 分 别 为 (x4,y4 ) 和 


(xs,ys)， 变 换 关 系 为 


xy =Ax+(1+m)x, cosO-(1+m)y, sind 


À (5) 
y, =Ay+(1+m)x, sin +(1+m) y, cos8 
转换 为 矩阵 形式 ， 
Ax 
x] [1 0 x, —y, || Ay| a=(1+m)cos@ kes 
y], [0 1 y, x Ja 'b=(1+m)sinð 6 
b 
误差 方程 为 
Ax 
V. |. 0 x ju Ay 2d = 
LA 0 1 y x ja LYs T Ya 
b 


Ht 


x Ax | 
其 中 ， M 为 转换 前 坐标 ; M 为 转换 后 坐标 ; p 为 平移 参数 ，(1+m) 为 缩放 参数 ; 


0 为 旋转 参数 。 理 论 上 ， 只 要 存在 两 组 转换 控制 点 ， 利 用 上 式 方 程 ， 可 联 立 解 出 一 组 转换 
四 参数 ;存在 3 个 或 3 个 以 上 的 控制 点 时 , 可 利用 最 小 二 乘法 求解 四 参数 X = (Ax, Ay m, 0)" ; 
然后 将 待 匹配 图 像 转换 到 参考 图 像 。 

需要 注意 ,在 视 场 平面 坐标 转换 时 ， 需 要 对 误 配 点 对 进行 筛选 ， 通 过 计算 控制 点 间 的 距 
离 。 在 实际 匹配 中 ， 如 果 控 制 点 之 间 的 距离 大 于 5 像 元 ， 即 2. 5” 时 ， 则 认为 该 点 为 误 配 点 
对 ， 将 其 吻 除 ， 根 据 入 选 后 的 控制 点 求解 的 RST 参数 ， 重 新 对 SDO/HMI 和 IM 太阳 望远镜 像 
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进行 RST 迭代 ， 计 算 RST 参数 ， 最 后 得 到 迭代 完成 后 的 RST 转换 参数 。 
3 实验 结果 与 分 析 


本 文 实验 中 ， 分 块 尺 寸 设置 为 30 像 元 *30 像 元 ， 块 间隔 为 5 像 元 。 图 3 为 利用 上 述 方 
法 ， 将 1M 太阳 望远镜 高 分 辨 像 和 SDO/HMI 全 日 面 像 配 准 后 重 关 区 域 的 对 比 图 ， | 


ee : 


a IM 太阳 望远镜 像 b SDO/HMI 对 应 的 区 域 
图 3 IM 太阳 望远镜 高 分 辨 像 和 SDO/HMI 全 日 面 像 配 准 后 重 闪 区 域 的 对 比 图 


Fig.3 Comparison overlapping area of 1M 太阳 望远镜 high-resolution image and full-disk Solar image after 
registration 
从 图 3 可 以 看 到 , 配 准 后 图 像 的 对 应 区 域 除 清晰 度 明 显 不 同 外 , 其 他 特征 已 经 完全 匹配 。 
为 了 更 好 地 观察 配 准 效果 ， 将 配 准 区 域 琶 加 后 ， 对 边界 区 域 进 行 放 大 ， 如 图 4， 可 以 看 到 边 
界 细节 的 信息 匹配 度 良好 ， 和 气孔 、 米 粒 、 米 粒 瞳 径 均 较 好 的 衔接 。 
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图 4 配 准 区 域 的 边界 细节 


Fig.4 The edge details of the registration area 


根据 配 准 结果 和 SDO/HMI 的 标准 参数 ， 表 1 显示 了 推算 的 这 张 IM 太阳 望远镜 图 像 PO 
角 和 像 元 比例 尺 ， 以 及 视 场 中 心 在 SDO/HNI 全 日 面 像 中 的 位 置 。 
表 1 配 准 结果 


Tab. 1 Registration results 


Zl] 从 
Pomme 比例 尺 some TAEAE Qus REIR 
= B 标 匹配 点 个 数 2 
终 配 ) 结 ( 度 ) (arcsecond 像 元 ) ( 像 元 ) (ix) (SDO 像 元 ) 
果 
-12. 95° 0. 04223 1524. 163 2686. 251 775 0.31 像 元 


从 表 1 可 以 知道 平均 残 差 为 0. 31SD0/HMI 像 元 ， 意 味 着 拟 合 结果 的 偏差 约 为 0. 15 ”。 
5 显示 了 这 些 残 差 的 空间 分 布 。 从 图 5 可 以 看 出 ， 残 差分 布 并 不 是 完全 随机 ， 显 示 了 1M 
太阳 望远镜 像 存 在 畸变 。 这 是 由 于 地 面 高 分 辨 观测 图 像 重建 后 , 依然 还 有 一 些 低 频 的 大 气 拌 
动 残余 未 能 消除 ,而 这 种 原因 限制 了 最 终 拟 合 结果 的 精度 。 同 时 由 于 两 个 图 像 的 观测 时 间 并 
不 完全 同步 ， 有 6 s 的 间隔 ， 因 此 有 些 太 阳 图 像 特征 已 经 发 生 了 改变 。 


Der e 


h 


WM cU T 


Mov e. 4 


AX LEM 


图 5 配 准 后 的 向 量 场 


Fig.5 The vector field after registration 


对 多 张 JM 太 阳 望 远 镜 和 新 太阳 望远镜 观测 的 高 分 辨 像 和 SDO/HMI 全 日 面 像 进行 了 匹配 ， 
结果 见 表 2。 
表 2 新 太阳 望远镜 和 1M 太阳 望远镜 的 TiO 图 像 与 SD0/HMI 全 日 面 像 的 配 准 结果 
Tab. 2 Registration result of TiO image(〈 新 太阳 望远镜 and 1M 太阳 望远镜 ) and SDO/HMI 
full-disk Solar image 


Ej SDO/HMI 的 FRET A AR 


FF "5 —_— TE — Ajs TI SH ` M 
$s 望远镜 高 分 辨 像 观测 时 间 ‘pies a 
观测 间隔 (SDO 像 元 ) 
1 新 太阳 望远镜 2014-08-05 22:29:59 ls 0. 33 
2 新 太阳 望远镜 2014-08-05 16:50:13 2 s 0. 43 
3 新 太阳 望远镜 2014-08-01 16:53:16 ls 0. 45 
4 新 太阳 望远镜 2012-07-05 16:36:28 Ts 0. 42 
1M 太阳 望 远 
5 M 2013-07-15 07:29:09 6 s 0.31 
镜 
1M 太阳 望 远 
6 H 2012-10-29 06:03:37 8 s 0. 50 
镜 
1M 太阳 望 远 
7 M 2013-06-12 07:45:34 lls 0. 26 
24 


从 上 述 结果 可 以 看 到 ， 最 大 的 匹配 残 差 约 0.5 像 元 ， 相 当 于 0.25" . 
实验 过 程 中 发 现 , 不 同 参 数 的 选取 会 对 算法 产生 一 定 的 影响 ， 具 体 见 图 6- 图 8, 


->- 30 degrees 
-0- 35 degrees | 
-*- 40 degrees 
-+ 45 degrees "| 


Residual / Pixel 


图 6 不 同 的 初始 参数 
Fig.6 The different initial parameters 
Al 6 显示 了 初步 估计 的 参数 对 迭代 次 数 的 影响 。 图 6 中 , SAW (绿色 )、 空 心 圆 型 CGE 
E)N KFE AEL HFE RE) 分 别 代表 估计 角度 为 30”、35”、40”、45” 时 的 
夺 代 过 程 。 可 以 看 出 ， 随 着 初 值 误差 增 大 ， 为 了 修正 误差 使 结果 收敛 在 合理 的 范围 内 ， 算 法 
的 迭代 次 数 相 应 增加 。 当 初始 估计 角度 达到 45” (此 时 的 误差 约 为 30”) IN, ARRA 
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6 次 增加 到 42 次 才能 收敛 。 


-+- 5 Pixel 


-0-10Pixel | 


图 7 不 同 的 控制 点 筛选 间距 
Fig.7 The different filter distance between control points 

图 7 显示 了 筛选 控制 点 间距 离 的 设置 对 算法 收敛 速度 的 影响 。 图 7 中 , 十 字 型 (红色 )、 

空心 加 型 (黑色 ) 分 别 代表 筛选 阔 值 为 5 Pixel 和 10 Pixel 的 迭代 过 程 。 可 以 看 出 ， 初 值 

误差 较 大 时 ， 高 阔 值 的 迭代 次 数 更 少 ， 意 味 着 算法 收敛 更 快 ， 反 之 低 冰 值 效果 更 理想 。 当 角 

度 误差 在 30” 附 近 时 , PERERA 28 次 ， 高 阔 值 更 优 。 角 度 误 差 在 10” 附近 

时 ， 和 迭代 次 数 仅 相差 OK, REEE. 


-*- Sigma=0.5 
-0- Sigma=1.0 
-+- Sigma=1.5 


Residual / Pixel 
o 


15 2 
Iterations 


图 8 不 同 的 平滑 参数 
Fig.8 The different filtering parameters 
图 8 Shas SATA BOT A SERI 图 8 中 , KFE CSR), toO A ( 蓝 
色 )、 十 字 型 〈 红 色 ) 分 别 表示 平滑 参数 为 0.5、1、1.5 时 的 迭代 过 程 。 可 以 看 出 ， 平 滑 参 
数 增 大 使 迭代 次 数 增长 ， 并且 拟 合 结果 有 一 定 的 降低 。 当 平滑 参数 为 0. 5 时 ， 算 法 达到 收敛 
所 需 的 迭代 次 数 和 结果 较 优 于 其 它 情况 。 


4 结论 


本 文 提出 了 一 种 针对 太阳 高 分 辨 像 和 全 日 面 进 行 视 场 匹 配 的 算法 , 从 而 解决 了 不 同 空间 
分 辨 率 下 、 上 共有 旋转 角 的 太阳 观测 图 像 匹 配 问 题 。 算法 的 基本 思想 是 结合 区 域 相关 匹配 和 控 
制 点 最 小 二 乘 求 解 的 方法 。1M 太阳 望远镜 和 新 太阳 望远镜 高 分 辨 图 像 与 SDO/HMI 全 日 面 像 
的 视 场 匹配 实验 都 证 明了 两 者 的 图 像 拟 合 结果 的 偏差 优 于 0.25”。 通 过 匹配 可 以 精确 测量 
高 分 辨 像 的 PO 角 ， 像 元 比例 尺 等 。 

这 套 方 法 要 计算 大 量 的 小 区 域 相 关 ， 同 时 进行 迭代 计算 ， 因 此 ， 时 间 开 销 较 大 ， 以 及 不 


同 参数 对 迭代 次 数 的 影响 ， 是 在 方法 应 用 时 需要 注意 的 。 除 此 之 外 ,小 尺寸 子 块 会 导致 匹配 
点 数 增加 ,计算 量 增 大 ， 拟 合 精 度 会 存在 一 定 的 下 降 ; 而 大 尺寸 虽然 对 精度 有 一 定 提升 , 但 
是 会 导致 匹配 区 域 变 小 , 从 而 会 影响 整 张 图 像 的 匹配 效果 。 而 子 块 间距 更 多 的 影响 计算 速度 。 
结合 实际 情况 ， 我 们 认为 子 块 尺寸 在 30 像 元 、 块 间距 5 像 元 附近 ， 算 法 的 稳定 性 较 合 理 。 
在 应 用 中 ， 可 根据 实际 情况 对 参数 进行 调整 。 
此 外 , 将 高 分 辨 Ti0 像 通过 与 SDO/HMI 全 日 面 连续 谱 图 像 的 高 精度 匹配 , 也 意味 着 可 以 
转换 为 和 任意 标准 的 全 日 面 像 进行 匹配 ， 如 SDO/HMI 磁 图 、SDO/AIA 各 波段 观测 ， 甚 至 和 色 
球 全 日 面 像 进行 匹配 。 这 意味 着 可 以 在 0.25 ”的 偏差 下 ， 研 究 小 尺度 结构 在 各 个 波段 ， 以 
及 在 磁场 上 的 对 应 情况 。 

这 套 方案 可 以 用 于 任何 具有 相似 图 像 结 构 的 太阳 图 像 匹 配 , 如 TiO Gband, Ha 偏 带 、 
10830; 也 可 以 用 于 Ha 线 心 图 像 和 GONG 全 日 面 像 的 配 准 。 因 此 具有 非常 多 的 应 用 前 景 。 
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High-accuracy Registration Method of Solar High-resolution 
Observation Images and Full-disk Solar Images 
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Abstract: In the research of the Solar observation, the registration between high-resolution 
observation images and full-disk Solar images is a greatly meaningful work. However, it is 
difficult to register with high accuracy due to the rotation, scaling and translation between them. 
This paper presents an image registration method combining local statistical information and 
control point matching: At first, between the high-resolution observation image and the full-disk 
Solar image is preliminary estimated about the orientation, scale, and positional parameters. Then 
the images are pretreated based on the estimated parameters. And then the pretreatment of the field 
of view is divided into a large number of overlapping local regions equally. Each corresponding 
local area is determined by the Correlation matching on the full-disk after preprocessing. Next, the 
sub-pixel offset between each pair of local regions is then measured, and the coordinate positions 
of the feature point are determined according to the sub-pixel offset between each pair of local 
regions, which is used as the feature control point in the point matching. Finally, the conversion 
equations of the affine transformation are established according to the control point, then use the 
least squares to solve the entire field of view of the transformation parameters. The image is 
reiterated according to the parameters to be solved, and the transformation parameters are obtained 
after the final iteration to register images. By registering the high-resolution observation images 
and the full-disk Solar SDO/HMI continuous spectrum, the deviation of fitting results is within 
0.25 arc-seconds. 


Key words: Solar Image Registration; High-resolution Image; Full-disk Solar Images; 


